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Ultracold fermions in a one-dimensional, spin-dependent, optical lattice are described by a non¬ 
standard Hubbard model with next-nearest-neighbour correlated hopping. Periodic driving of the 
lattice allows wide tuning of the system parameters. We solve this model exactly for a special value 
of the correlated hopping. The solution reveals the general properties of this system for arbitrary 
filling: exact and asymmetric spin-charge separation, a gapless spectrum of lowest energy excitations, 
a spin-gap, which may be interpreted in terms of collective hole pairing and a non-vanishing den 
Nijs-Rommelse type string correlator. Numerical studies away from the integrable point show the 
persistence of both long range string order and spin-gap. 


Introduction - The identification and description 
of interesting collective states of many particles is the 
essence of contemporary many body physics. The study 
of one-dimensional (1-d) systems holds a special place in 
this effort [T], where effects of enhanced quantum fluc¬ 
tuations are particularly important. These systems are 
amenable to various exact or quasi-exact analytical and 
precise numerical methods which reveal surprises and ex¬ 
otic physics, e.g. Luttinger liquids and associated spin- 
charge separation ma, fractional excitations [5], non¬ 
local symmetry breaking captured by string order pa¬ 
rameters [71IH] , etc. Due to recent progress in ultra-cold 
quantum gas experiments, simulations of various theo¬ 
retical models with unprecedented control are becoming 
possible. In particular, non-standard extensions of the 
paradigmatic Hubbard models are highly desired |5]- A 
notable example concerns correlated hopping (CH) con¬ 
sidered already in the context of multi-orbital physics 
[lOj . dipolar interactions m and time-periodic modula¬ 
tion of on-site interactions [12j . 

In this Letter, we investigate a 1-d fermionic model 
with an extended CH spanning next-nearest nieghbours. 
This model precisely describes a system of ultracold 
fermions loaded onto an optical lattice made of two 
shifted, spin dependent sublattices. This structure makes 
the CH naturally stronger and/or comparable to normal 
hopping. A suitable periodic tilting of the lattice, effec¬ 
tively weakens the interactions between fermions and al¬ 
lows to access the dynamics due to the strong CH. We de¬ 
scribe the ensuing physics analytically for a special value 
of the CH and numerically in the general case. We pre¬ 
dict numerous exotic properties of the system: exact and 
asymmetric spin-charge separation, a gapless spectrum of 
lowest energy excitations, a spin-gap associated with col¬ 
lective hole pairing, non-vanishing string correlator and 
transition to a simple ferromagnetic state. 

The standard, nearest-neighbour CH originates from 
inter-particle interactions, as does the diagonal on-site 
interaction [13] . Such CH was intensively studied for 
strongly correlated electrons, in particular, as a mecha¬ 


nism for unconventional (non-phonon mediated) super¬ 
conductivity mus], especially at weak to moderate val¬ 
ues of the CH integral X. Large values of X comparable 
to normal hopping J block certain tunnelling processes 
resembling effects of large repulsive Hubbard interaction 
[Mill]- For extreme correlations X = J, fermion motion 
is kinetically constrained so that the on-site interaction 
term U commutes with the rest of the Hubbard Hamilto¬ 
nian. This additional symmetry facilitates an exact so¬ 
lution in 1-d nBHii], revealing several non-perturbative 
effects such as a paramagnetic Mott transition (at half¬ 
filling) and exclusion statistics jH]. Furthermore, an in¬ 
commensurate singlet superconducting phase appears in 
the proximity of the X = J limit |23l |2l] . Separation of 
lattices available for hopping of the two species of spins 
leads to a model with a different kind of CH, extended 
over three-sites. We study this model in this Letter re¬ 
vealing markedly different phenomenology as compared 
to the standard CH Hubbard model. 

The model - Consider an optical lattice setup with a 
red and blue detuned sublattice offset by half a lattice 
constant (for related experiments on spin-dependent lat¬ 
tices, see |75H77] 1. each loaded with only a single fermion 
spin species t or /. Equivalently, one could load two 
species of Tonks-Cirardeau bosons [28l - l^ . Each spin 
species only hops on its own sublattice and different 
fermions undergo contact interactions. As there is no 
mixing between sublattices, the spin index can be sup¬ 
pressed yielding the Hamiltonian: 

H = + Ani+i)(cJcj +2 + h.c) + ^ VuiUi+i, 

i i 

( 1 ) 

where odd (even) sites correspond to the blue (red) sub¬ 
lattice. The hopping J > 0 connects next nearest neigh¬ 
bours. The leading term stemming from the contact in¬ 
teraction is the nearest neighbour repulsion V 15> J, fol¬ 
lowed by the CH to a next nearest neighbour > 0 con¬ 
ditioned on the intervening site occupancy (Figj^ details 
in Supplemental Materials - SM). Interestingly, although 
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FIG. 1. (Color online) Model parameters vs. lattice depths 
Vo for a scattering length = 300ao (Bohr radii) and har¬ 
monic trap energy hutm-p = {Eh - recoil energy). For 

shallow lattices, X « 1.6J and decreases to X = J for a deep 
Vo « 14.25i?H attaining a value X ~ 0.88J for Vo = 18Er. 


FIG. 2. (Color online) Mapping to spinless fermions in the 
kinetically constrained model {X = J,V = 0) when different 
spins cannot pass each other. Occupied sites are shaded. Each 
pair of empty sites between different spin consecutive fermions 
is an empty site in the charge chain. Empty sites on a single 
chain between two closest same species fermions are empty 
sites in the charge chain. 


in usual scenarios CH terms are subleading in comparison 
to normal hopping [5], the separation of the sublattices 
yields a novel, viable path to a strong CH term X ~ J. 

Fast, periodic tilting of the optical lattice described by 
Hs{t) = Acosojt'^- jrij [211133], at frequency close to 
resonance with the interaction: V = 2a; -I- <5 with detun¬ 
ing 6, allows control over the full range of properties of 
the system. The long-time dynamics of fermions is de¬ 
scribed by the Hamiltonian H{t) = H + Hg{t) in the 
frame moving with the interaction and the lattice tilting: 
H' = UHW + iilW , with U = exp(z2a;t njUj+i + 

i Jq Hs{t')dt'). As shown in the SM, choosing the tilting 
amplitude A so that Jo(2A/cj) = Jii^Alio), where Jn{x) 
are the n-th order Bessel functions of the first kind, the ef¬ 
fective Hamiltonian H' is time-independent and takes ex¬ 
actly the form of the original CH model ([^ with dressed 
parameters J —>■ JJo{2A/uj),X —>• XJq{2A/uj),V —>■ 6. 
So, although the bare interaction is large for fixed ratio 
X/J, the driving makes relevant the complete range of 
parameters of the model Q from V = 0 (exact resonance 
at 5 = 0) to J/V = 0, controlled via the detuning and the 
dressing parameter Jo{2A/oj). We elucidate the model’s 
properties below. 

Exact solution at X = J,V = 0 - At this point, the 
model has an important symmetry, as the free and corre¬ 
lated hopping cancel for two nearest neighbour fermions. 
So, different fermions are forbidden to pass each other 
and the sequence of f and fermions is conserved under 
open boundary conditions (OBC). The Hamiltonian can 
be diagonalized for each fermion sequence separately and, 
as shown below, fermions may be completely stripped of 
their spin. The fixed sequence in the original lattice con¬ 
strains the space available for the motion of fermions. 
Thus while the original system consists of two coupled 
chains of length L, the spinless fermion degrees of free¬ 
dom live on a charge-chain of length Lgff i which remark¬ 
ably depends on the fermion sequence. 

We restrict the total number of fermions Nt < L, so 
as to allow the realization of all possible spin sequences. 
Note that the distance (number of empty sites) on the 
original lattice d between two consecutive fermions is al¬ 
ways even (odd) if they belong to different (same) species. 


with the minimal distance 0 (1) respectively (see Fig. [^. 
The distance of the first (last) fermion from the leftmost 
(rightmost) edge of the lattice is always even for I (f) or 
odd for t (1) fermion respectively. The conservation of 
spin order leads to the conservation of order of parities 
of consecutive distances. For example (see Figj^ for 6 
fermions and the following spin order itiitt there are 
7 distances and their parity order is EEEOEOE. So, to 
specify any configuration, for a given spin order, one may 
use a set of effective distances des, where des = d/2 or 
des = (d — l)/2 for even or odd d respectively. The d^s 
are the number of empty sites in the charge lattice. Note 
that the sum of these effective distances Dgs is given by 
Defi = {Dq — Mo)/2 where Dq = 2L — Nt is the sum of 
the distances on the real lattice and Mq is the number of 
odd distances. This means that Dgg is a constant for a 
given spin order, and there is a 1-to-l mapping between 
configurations on the real lattice (for that order) and con¬ 
figurations of Nt spinless fermions on the effective lattice 
of length Leff = Des+Nt = L+{Nt—Mo)/2. Futhermore 
note that the number S = {Nt — Mo)/2 is the number of 
switches, i. e. the number of times a fermion is followed 
by a t fermion. Hence L < L^e = L-\-S < L + N^, where 
N^ is the number of fermions on the upper chain. 

The action of H (with H = 0) in the original lattice is 
equivalent to the action of nearest neighbour free-fermion 
hopping in the charge lattice, which immediately allows 
the identification of all eigenstates and energies of the 
system. In particular the lowest energy in a given spin 
order sector, in the thermodynamic limit, is 


Eo/L 


2J -Leff . 7T2pL 

-sin- 

TT L -Leff 


( 2 ) 


where p = Nt/{2L) is the original fermion density. The 
most interesting features of the solution are listed below: 

(i) The decoupling of the spin order from the charge 
degrees of freedom is an intricate manifestation of spin- 
charge separation ubiquitous in 1-dimensional settings 
[T]. The novelty here is that it is asymmetric: while 
charge excitations can be created without affecting the 
spin configuration, a change in the spin order explicitly 
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affects (via ) the charge energy spectrum. 

(ii) We restrict to an even number of fermions Nt = 
2N. The global ground state corresponds to spin order 
of the form | itit ■ • • it) • This sequence maximizes 
Leff for a given L and 2N as it corresponds to the low¬ 
est possible particle density on the charge lattice. Hence 
the spin part of the ground state wave function has per¬ 
fect antiferromagnetic order. Note that this order is not 
degenerate with respect to flipping all spins - this is a 
manifestation of the lack of reflection symmetry of the 
zigzag chain about a horizontal line through the middle 
of the chain. 

(iii) Antiferromagnetic order in the spin sector implies 
the existence of hidden non-local string order in the full 
lattice of moving particles |34j . where each fermion is 
followed by an f, with an arbitrary number of holes be¬ 
tween each pair. This is an analog of string order in 
spin-1 Heisenberg chains [7]. We stress that there is no 
local symmetry breaking of the density wave type accom¬ 
panying the formation of this string order - the ground 
state is translationally invariant (modulo OBC). Unlike 
for spin -1 chains, the low lying excitations are gapless (in 
the thermodynamic limit) corresponding to charge exci¬ 
tations above the Fermi sea at fixed spin order. 

(iv) Although the spectrum is gapless in each order 
sector there is a spin gap, i.e. the energy cost asociated 
with flipping one spin which shortens the charge lattice 
size by one site L^g = L + N — 1. All single particle 
energies are shifted due to the change in Leg leading, 
in the thermodynamic limit, to a non-zero spin gap for 
arbitrary Ailing densities. 
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FIG. 3. (Color online) All nearest neighbour holes are 
uniquely paired in the ideally string ordered configurations of 
fermions. So, the motion of fermions in the real lattice can be 
equivalently described as hole-pair hopping. (Below) Hole- 
pairs are single particles that hop on a dual lattice defined 
on the edges of the zig-zag chain (linear chain of squares). 
Shaded squares are bound hole-pairs. Uniqueness of pairing 
of holes (above) implies exclusion of simultaneous occupation 
of nearest neighbour sites on the dual lattice, i.e. infinite 
repulsion of paired holes on these sites. 


symmetry breaking but rather algebraically decaying su¬ 
perconducting correlations (SC). Explicitly, notice that 
the motion of paired holes in configurations as in Figj^ 
can be described by a spin-1/2 XXZ model at fixed mag¬ 
netization m = p/ 2 , defined on the edges of the zigzag 
chain: H = J2^ >/(cr|tTj+i -b h.c) + + l)^^-! + 1 ), 

where = 1 (— 1 ) corresponds to the occupation (lack) 
of a paired-hole and A = oo. The ground state is a Lut- 
tinger liquid and the SC in the real lattice correspond to 
the transverse spin correlation function in the dual one. 
The latter can be obtained using bosonization techniques 
(see Ch. 6 and Eq. (6.47) of Ref. [1]): 


As = Eo(Leff = L + N-l)- EoiLeg = L + N) 

2,7((p + l)si„(|m)-2xpcos(|™)) 

For X = J switching the order of a pair of nearest neigh¬ 
bour fermions in the ground state entails the same energy 
cost. There are N such degenerate spin order configura¬ 
tions. Interestingly, the alternative perfect antiferromag¬ 
netic configuration | titi • • ■ ti) obtained by flipping 
all spins in the ground state is also associated with a 
{L + N — l)-site charge lattice and so is exactly degen¬ 
erate with the lowest lying spin excitation. The As 7 ^ 0 
suggests a protecting mechanism for the string ordered 
phase in our system when moving away from the inte- 
grable case X = J, V = 0. 

(v) A revealing aspect of the ground state is that all 
holes are unambigiuosly bound into nearest neighbour 
pairs (Figj^. The particle motion corresponds to the 
coherent motion of paired-holes through the lattice (in 
the opposite direction to a particle). The finite spin gap 
As for all densities, which in the paired hole picture is 
the energy of breaking paired holes, suggests the ground 
state can be interpreted as an exotic hole superconductor. 
Since, the system is 1 -dimensional, there is no explicit 


(Ca;+lCa;CQC-|^)i.eal — {S (x)*S (0))dual 

.ls2K+l/(2K) /l^l/(2if) 

= C cos{2'Kmx) y-j -b C"cos(7ra:)(^—j , 

where the Luttinger parameter K increases from 1/4 for 
p —)• 0 to 1 for p —)• 1 [3T] . 

Finally, note that while the configurational mapping 
to the effective chain holds for A = J, U ^ 0, the solely 
interspecies interaction negates the mapping to a spinless 
fermion Hamiltonian in general. However, two spin or¬ 
ders are exceptions: the perfect antiferromagnetic order 
maps onto a. t — V spinless fermion model which as seen 
below is stable upto finite V, and the fully polarized spin 
state which stabilizes at large interactions. In particular, 
for U = 00 {J/V = 0), the largest possible obtainable 
length is Leg = L yielding a ferromagnetic (FM) phase 
separated ground state for equal number of f and ], spins: 

I U • • • itt • •. t). 

Robustness of string ordered phase - Intuitively, to 
induce a change in the described ideal string configura¬ 
tion for X = J, the energy penalty associated with V 0 
should be large enough so as to offset the kinetic energy 
gained via the introduction of defects in the spin order. 
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FIG. 4. (Color online) Exact Diagonalization (ED) (a,b,e,f) and TEBD algorithm |36i I37 | (c,d) results, (a) O vs F for 
X = J where 2N = 8,2L = 20 (black), 2N = 8, 2L = 24 (red),2A'^ = 6, 2L = 24 (blue), (b) Extrapolated spin gap to 
thermodynamic limit for p =1/2 (red) and p = 1/3 (blue) for V = 0. (c,d) Finite size scaling Os vs distance x under OBC for 
different densities and V = 0. Data are for X = l.OJ (diamonds), 0.6J (circles), 0.3J, (squares), O.OJ (triangles). Compared 
are lattice lengths 2L = 24 (blue), 2L = 48 (red) and 2L = 60 (green) for every X. (e-f) Os for on a lattice of length 2L = 20 
for OBC (e) and periodic boundary conditions (PBC) (f) for V = 0 and X = J, X = 0.5J, X = OJ from top to bottom. 
Oscillations in data in (c-e) are due to Friedel oscillations and OBC. These are diminished under PBC (f). 


In order to analyse this, we introduce the parameter: 

2L 

where Qjj- = nj'>z>j(l “ ”0 for / > j + 1 or Qjj, = 1 
where j,j' are subsequent sites. With an equal number 
of t and i fermions, O = 2N — 1 is maximal for ideal 
string ordered configurations, and minimal for the phase 
separated state with 0 = 3 — 2N. Ground state numeri¬ 
cal results (see FigQa)) show that for arbitrary densities, 
a large region of stability of ideal string order (and thus 
collective hole pairing) is followed by a first order tran¬ 
sition to the FM state for V < 2J. In some cases, the 
jump - in these finite size systems - is via an intermediate 
spin sequence. It is interesting that the simple variant of 
itinerant ferromagnetism unambiguously exists here for 
realistic parameters. In particular, without lattice driv¬ 
ing and fine tuning, the bare interaction (Fi g§ places 
the system in the FM phase. 

For X ^ J, the spin sequence is subject to dynamical 
changes. We restrict to X/J < 1 due to the symmetry 
X 2 J — X of the model Q [3T]. Without CH {X = 0, 
V yf: 0) an interaction driven, direct liquid to phase- 
separated jump transition for arbitrary fillings already 
exists [35]. So, we shall here not consider the similar ex¬ 
pected effect of F 7 ^ 0 but rather the fate of the string or¬ 
dered state for F = 0. First, note that the spin-gap As is 
always open for X > 0 (Fig. |^b)) indicating hole pairing 
and long range string correlations for arbitrary low val¬ 
ues of CH. To probe long range non-local hidden order for 
general X ^ J, we introduce a den Nijs-Rommelse type 
string parameter lilMj, Os- Grouping pairs of nearest 
neigbour sites into dimers (1, 2), (3,4),... (2L—1, 2L), de¬ 


fine the operator Sf = (n 2 i—n 2 i-i) with values +1,0, —1 
in analogy with a spin-1 Sz operator, so that 

Os(m,m') = (5'^exp (^fTT ^ (5) 

m<l<m' 

Figs. Qc-f), clearly show that Os attains a finite value 
for the ideally string ordered ground state at X = J. Re¬ 
ducing X leads to a diminishing, yet finite large distance 
saturation value of Os for X 0 with no evidence of 
destruction of long range string order before X = 0. We 
have checked that this holds for various densities p. 

Final comments - Long range string order is often 
associated with a symmetry protected topological phase 
characterized by existence of edge modes and topological 
degeneracy. Our model, under OBC, breaks inversion 
symmetry which has been shown recently to protect the 
existence of such features [35]. Indeed, the string ordered 
ground state is unique in the thermodynamic limit, since 
the spin gap separates it from the configuration with all 
spins flipped. This feature is distinct from the standard 
Haldane phase - there are no zero energy edge modes 
yielding ground state degeneracy. 

While string order parameters were first shown to char¬ 
acterize gapped spin-1 chains laiHj, recent studies have 
shown that specific non-local order parameters naturally 
describe a plethora of gapped phases without local or¬ 
der in 1-d systems. However, these commonly take the 
form of parity correlators that capture the confinement 
of quantum fluctuations e.g. in Mott insulator phases of 
bosons [39] and fermions as well as the Luther-Emery 
phases of the attractive Hubbard model [40l [41]. In- 
triguingly, however, critical string order (algebraically de¬ 
caying) was shown to underlie spin-charge separation in 
1-d systems surviving upto the limit of non-interacting 
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fermions [34) . It is interesting that the asymmetric type 
of spin-charge separation encountered in the model con¬ 
sidered in this Letter leads to a truly long range string 
order. Moreover, the lowest lying charge excitations in 
the system are gapless and our model thus constitutes a 
rare example of finite string order in a gapless phase (see 
also Refs. 0133]) of itinerant particles. 
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SUPPLEMENTAL MATERIALS: 


EFFECTIVE PARAMETERS OF THE 
CORRELATED HOPPING MODEL 


Our system consists of two-species of ultracold atoms 
trapped in a 1-d optical lattice. For one species, we con¬ 


sider a blue detuned laser and for the other a red-detuned 
laser to create the optical lattice. As a result, the respec¬ 
tive lattices for the species are shifted in space by a half¬ 
lattice constant a/2. The lattice potentials are therefore 
^iatt(^) = V^cos^ -Kx/a and = 1^ cos^ (7rx/a-I-(/) 

where the red and blue-detuned case corresponds to 
(j) = t: 12. The corresponding Hamiltonians reads. 



■ (f 

2ma dx"^ 





(SI) 


where 'hj.(x),are the creation and annihila¬ 
tion field operators for the species cr The 

first term in [..] describes the single particle kinetic 
and potential energy while the second term describes 
the one-dimensional contact interaction with strength 
5iD = reduced mass m, = 

2m^mi/{m^ + mi), scattering length Og and the trans¬ 
verse oscillator length Otrap and constant C « 1.46 


[T]. Next we expand the field operators, = 

a;i^(a;)cr|_^, where v denotes the sites for the species 
a =t, i, (x) are the respective Wannier functions and 
aj are the creation operator for the species a at site ia- 
The contact interaction leads to both nearest neighbour 
repulsion and correlated hopping in the Wannier basis. 
Indeed, the corresponding tight-binding Hamltonian in 
terms of the local operators reads: 






+1 -I- (Tj -I- h.c) + X + 






i-j- +1 


)nl 


(S2) 


To express the parameters, first we assume that = 
mi = m and the lattice depths are same, Vi = Vi = 
Vo- Moreover, we express the lengths x ^ xja and the 
energy with respect to recoil energy Er = I2ma^. 
This implies that Ji = Ji = J > 0 and for Wannier 
functions we use the convention, ujiy.{x) = uj{x — Xi) and 
ojiy^{x) = oj{x — Xi — 1/2). With this convention Eq.(S2) 
is equivalent to the concisely written model Eq.(l) in the 
main text with suppressed spin indices. The parameters 
are now expressed as, 



^ trapping potential in 

Fig.l of the main text corresponds to atrap/a « 0.14. 


PERIODIC LATTICE MODULATION AND 
EFFECTIVE TIME INDEPENDENT 
HAMILTONIAN 

We consider the system described in the previous sec¬ 
tion to be subject to shaking of the optical lattices with 
driving frequency ui and tilting amplitude A. The time- 
dependent Hamiltonian is then given by: 

H{t) = H + Acos{ujt) 

3 

where we represent the Hamiltonians in concise form cor¬ 
responding to Eq.(l), i.e. j enumerates consecutive sites 
on the zigzag chain. 

The interaction strength V is the largest energy scale 
in the system. We choose the driving frequency to sat¬ 
isfy the quasi-resonance condition V = 2u) -\- 5, where 
5 is the (small) detuning. In order to eliminate the 
fast degrees of freedom in the problem, we transform 
to the co-rotating frame via the unitary transformation 
U = e'yii)(i2iut^^rLjnj+i+iA/uj juj) to obtain 
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the effective Hamiltonian H' = UHW + itlW: 


time dependent function 


iJ' = ^ Qi(t){-J + Xni)c\ci +2 + h.c. + (5^ rurii+i 

i i 

(S5) 

where the hopping and CH terms are modulated by the 


Giit) = exp 



2A 

i — smwt 


(S6) 


The occupation number operator is a projector = n, 
so ex\)(iBtn) = (1 — n) + nex:p(iBt), and therefore 


exp(—— 712 )) = (1 — ni — 77-2 + 2nin2) + ni exp(—tiJt) + 77,2 ejq){iBt) — nin 2 {exp{iBt) + exp(—tHt)). (S7) 


Using this and the Jacobi-Anger identity, 

UJ 


( \ /2v4\ 

exp —sina;tj = 2_^ —j exp(zmu;i) (S8) 


m— — c>Q 


we obtain the following time independent effective mod¬ 
ulation G on averaging out the fast oscillating terms: 


G^ 


— (1 — 77i+3 ~ Tli-i + 2r7i+377i_ 



+{nt+3 + rii-i - 


277 ^+ 377^-1 



(S9) 


after applying the relation J 2 {x) = J- 2 {x) for the Bessel 
functions of the first kind. The time averaging procedure 
can be carried out as long as w J,X,6 which is the case 
here (see Fig.l in main text). In general, the lattice mod¬ 
ulation generates CH terms, through Gi, beyond those 
considered in the paper. However, choosing the tilting 
amplitude A* such that Jo(2A*/a;) = J2{2A*/co), elim¬ 
inates these completely and one obtains the CH model 
Eq.(l) with the bare hopping and CH terms renormalized 
by ^ = Jo{2A*/uj), i.e. 



J -\- Xni)c\ci +2 + h.c. -I- 5 ^ riirii+i 

* (SIO) 


The hopping renormalization is given by Jq( 2A*juj ~ 
1.84118) = 0.316 at the first intersection point of the 
Bessel functions. The residual nearest neighbour interac¬ 
tion now stems solely from the detuning b and in partic¬ 
ular can be chosen to be equal to zero as well as non-zero 
values - both cases are considered in the main text. No¬ 
tice that choosing larger amplitudes A* corresponding to 
subsequent intersection points leads to a fast diminishing 
of the effective normal and correlated hoppings thus al¬ 
lowing to access smaller effective hopping to interaction 
ratios. 
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